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Abstract 

We compute the quark-antiquark potential employing the continuum action for QCD-like random 
lattice strings proposed by Siegel. The model leads to a potential similar to those obtained from 
Nambu-Goto string theory but has some modifications which we interpret as velocity dependent 
contributions. We also propose to add extra terms in the action which lead to physically interesting 
propagators for partons for the infrared region. 
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I. INTRODUCTION 

Despite the currently established status of QCD as the theory for strong interactions, 
due to its non-perturbative nature in the infrared regime, it is desirable to seek a simpler 
description of low-energy hadron spectrum in terms of phenomenological models. Often 
these models involve some type of non-trivial classical structure like bags, strings or in- 
stantons. Amogn these strings were historically introduced to explain the so-called Regge 
behavior and duality symmetry demonstrated by the hadron-hadron scattering amplitudes 
at high energies. However, quantization of the strings as fundamental hadronic objects 
suffered from several severe drawbacks: the presence of unphysical particles in the spectrum 
( tachyons, massless spin-2 particles etc. ) as well as requiring a 2 6- dimensional spacetime. 
The subsequent advent of the QCD has since then pushed the hadronic string model into 
its current status as a phenomenological attempt to explain the hadronic scattering. Fun- 
damental strings have instead taken up the role of the prime model for the unification of 
fundamental interactions. 

Despite the false start, string theory received some theoretical justification from QCD 
when 't Hooft [|IJ showed that SU (N) gauge theories behaved as an effective planar ( hence 
a two-dimensional ) theory in the limit iV — > oo, provided color is confined. An earlier 
connection between string theory and large gauge theories was provided by the Nambu 
conjecture relating the Wilson loop operator ( defined in gauge theoretic context) defined 
over the closed loop C to the string theory defined on the world sheet £ bounded C ( cf. 



Here the string is assumed to be moving freely in a flat spacetime. In the recent years, this 
correspondence has been put on firmer foundations by Maldacena|2[ and Witten || which 
currently goes by the name of AdS/CFT correspondence M where the low-energy string 
theory is represented by five dimensional AdS supergravity theory and the corresponding 
four- dimensional gauge theory lives on the minkowskian ( conformal ) boundary of the AdS 
spacetime. 

An alternative path, explored by Nielsen and others ||, is to show the amplitudes of 
the dual models can be reproduced from scalar field theories ( representing the "partons" 
) defined on two-dimensional spacetime using Gaussian propagators for the scalar fields. 
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FIG. 1: String worldsheet £ bounded by the loop C over which supports the Wilson loop 



Though such propagators provide the field theoretic reason behind the better high-energy 
behavior of the scattering amplitudes of strings, the underlying field theory needs to be 
a non-polynomial theory. Also, it is known that Gaussian propagators behave as massive 

_ p 2 

propagators for low value of momentum, A 2 e _ A 7 ~ A a +p 2 wmcn contradicts our knowledge 
of perturbative QCD. This tells us that such "partonic" models do not corresponds to QCD 
string. 

Therefore, it is natural to attempt to construct a string action action whose "partons" 
are massless. Indeed, sometime ago Siegel has derived such an action |J. This action is 
particularly interesting from the fact that it correctly predicts the spacetime dimension to 
be four by demanding T-duality invariance of the string action. However, to explain the 
asymptotic freedom one has to assume that the interactions in this theory be described by 
the wrong sign 4 theory. Though, this theory will suffer from the problems with vacuum 
instabilities we will not concern ourselves here. 

An important feature that emerges from string models is that a linearly confining potential 
is obtained in the leading order from string actions proportional to the worldsheet area : 

V{R) = tR + ■ ■ ■ (1.2) 

where r is the string tension and the sub-leading order terms are represented by the ellipsis. 

In this brief report, we attempt to check whether this theory manages to reproduce 
the typical potential between quark-antiquark calculated from string theories. Unlike the 
Nambu or Polyakov strings the Siegel string involves a traceless symmetric matrix apart 
from the regular worldsheet metric. We find that the off-diagonal elements of this matrix 



contributes to the one-loop effective potential. We suggest that these off-diagonal elements 
are related to the velocity of the endpoints of the strings where the quarks reside. However, 
the matrix elements are not fixed in this theory. This problem is quite similar to the runaway 
problem of the moduli fields in fundamental string theory. However, unlike the fundamental 
string theory writing a potential for this matrix field is not prohibited anything like general 
coordinate invariance or supersymmetry. We argue that a realistic potential modifies the 
partonic propagators in the infrared regime which is reasonable for a string picture. 

The paper is arranged as follows. We first sketch Siegel's construction of the action for 
QCD-like strings in section 2. In section 3, we carry out the one-loop computation of the 
quark- ant iquark potential using Siegel's action following the procedures of Liischer et. al. 
and Alvarez ||. We contrast our results with the well-known results JF] in the literature. 
Finally, we comment on our results and on its physical implications. 



II. SIEGEL'S ACTION FOR QCD-LIKE STRINGS 

Siegel's construction of the action stems from the old observation due to Nielsen and 
others || that one can reproduce the amplitudes of the dual models by working with a two 
dimensional scalar field theory whose propagators are Gaussian in the momentum space. 
Therefore, it is interesting to search for a string theory which "originates" from a theory 
of massless propagators. Following the random triangulation approach for low-dimensional 
strings, we consider the worldsheet to be represented by a random lattice. The continuum 
limit will be recovered when the number of points on the lattice N — > oo. 

One begins with the Fock-Schwinger-DeWitt propertime parameterization of the massless 
propagators, namely 

i r°° 
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Consider a Feynman graph involving massless scalar fields on a random lattice ( cf. Fig |2] ) 
whose points are labelled by integers % — 1, 2, • • • , N with N — > oo. The propagators denoted 
by pij , being two point functions, are elements of a N x iV matrix and supported on the 
link connecting the lattice points i and j . The amplitude in momentum space associated 
with this graph will be given by 

A~ f U dD P^U S CLP^ (2-2) 

•> <ij> ' v i j 



FIG. 2: A random lattice graph which makes up the string worldsheet 
where the delta function simply enforce the momentum conservation at all vertices. Using 



the Schwinger parameterization (|2.1|) in (|2.2|) we rewrite the amplitude as 



a ~ I n <i "i > u iir u ( ; I "n dx ^ EjpiiXi ( 2 - 3 ) 

<ij> i 

where we have used a matrix-valued Schwinger parameter. The momentum conservation 
equation at the vertices can be solved by going to the loop representation for the graph, i.e. 
by writing 

Pij = Pu = kj-kj (2.4) 

where kj, kj are the loop momenta running over the ith and Jth loop, whose intersection 
is the link represented by the pair ij. momenta to obtain 

A ~ [ Y[dTjje- TIj{kl - kj)2 (2.5) 
J i 

where we have switched to the loop representation. 

To obtain the continuum limit of this action one has to note that there are only two 
independent directions on the string worldsheet one demands that the matrix has two degrees 
of freedom - so this matrix must be a traceless symmetric matrix. 

The final result is the Siegel's action is of the form 



S 



J d 2 a [r mn d m X ■ d n X + Xg mn r mn + • • •] (2.6) 



where g mn is the worldsheet metric used for taking the trace of r and the last term imposes 
the tracelessness condition via the Lagrange multiplier A. 
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FIG. 3: Quark-Antiquark Loop used for calculating the Quark-antiquark Potential 
III. CALCULATION OF THE QUARK-ANTIQUARK POTENTIAL 

The calculation of the quark-antiquark potential using the string models was pioneered 
by Liischer, Symanzik and Weisz J7| using the Nambu-Goto action. Further improvements 
of this result were done by Alvarez || and Arvis ||. 

The calculation proceeds via considering a quark-antiquark virtual rectangular loop where 
the separation between the quark-antiquark pair is R and the pair propagates for time T as 
in figure |3|. In the confining phase, the Wilson criterion gives 

< W{C) >~ e- TV ^ (3.1) 

Using the Nambu conjecture one can relate the string action with the quark-antiquark 
potential: 

V(R) = ±S(£) (3.2) 

Thus in flat spacetimes when the string action is proportional to the worldsheet ( hence 
the loop) area one recovers a linear potential in the tree level. However, at short distances 
the quantum fluctuations can have important consequences for instance determining the 



deconfmement behavior [10]. Therefore it is important to calculate the one-loop contribution 



to the quark-antiquark potential. 



We begin our calculations by embedding the worldsheet in the spacetime as 

T = t = X°, a = X l 



(3.3) 




The rest of the coordinates act as transverse excitations on the worldsheet and contribute 
to the effective potential via quantum fluctuations. 
Let us parameterize the matrix r mn as follows 

(3.4) 

which is traceless with respect to a Minkowskian metric. As briefly mentioned earlier, the 
expectation values of the fields a, b are not fixed due to the absence of any classical potential 
for them. However, following the cr-model approach to string theory we assign a constant 
vacuum expectation value for a, b ( this type of approximation was also used by Alvarez || 

)• 

Using the embedding ( |3.3| ) the action now reads 

S = 2aTR + [ d 2 aT mn d m (j) ■d n <f>+--- (3.5) 



One can see that the diagonal matrix element a is gives rise to the tree level string tension. 
However, the transverse fluctuations also contribute to the string potential which is to be 
interpreted as the contribution from the zero point fluctuations about the classical configu- 
ration (|3.3| ) ( Casimir energy ). If we integrate out the transverse fields the "effective action" 
reads 

S = 2aTR + (jP ~ 2) Tr In [d m r mn d n } (3.6) 

To find a real effective action we have to guarantee that the eigenvalues of the operator 
appearing in the second term of ( |3.6|) be real - hence the operator be self-adjoint. This can 



be done by choosing appropriate boundary conditions |L3[ as the worldsheet has boundaries. 



The self-adjointness condition for the operator appearing in (|3.6| ) would read as 



{ X ,V^) = {V X ,4>) V^,xeM (3.7) 



where M is the domain of the operator P[13|: 



V = d m {r mn d n ) (3.8) 



specified by the chosen boundary condition. In equation (|3.7Q the inner product is defined 
is the standard fashion, namely, 



(r],0= f R dar,\a)i{a) 
Jo 



(3.9) 



We are interested in the stationary configurations of the string and accordingly seek 
stationary solutions of the form 

i/j{a,t) = i){a)e~ iujt (3.10) 

Restricted to this type of subspace, the operator V takes the form 

= [a{dl - lo 2 ) - 2ibud a ] (3.11) 

Due to the presence of the worldsheet boundaries, the condition Q3.71 ) for the operator 
V w now translates into the vanishing of the following anti-hermitian bilinear form 



(3.12) 



where the vectors x, £ are of the form 



(3.13) 



and 



where tana 



CUx = o 

4>(R) 



2itana 1 
-10 
-2itana -1 
10 
- G ( — j , j ). The limits follow from the conditions that the eigenvalues of 



A 



(3.14) 



the matrix r must be positive and accordingly det r = (a 2 — b 2 ) > 0. Therefore - < 1. 

The eigenvalues \ii and the corresponding eigenvectors with % = 1 
A are as follows: 



• 4 of the matrix 
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(3.15) 



fj, 2 — i tan( 
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cos # + 



sin(f + I) 



(3.16) 
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(3.17) 
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(3.18) 



the eigenvectors have been normalized here. 

Domains over which the operator T>^ is self-adjoint can now be found by taking various 
combinations of the eigenvectors ( |3. 15 , 3. 16 , 3.17 . 3.18 ) with relative U(l) phase rotations 
between them. Therefore, any function^ belonging to the domain M will be given by a 
general linear combination of the form 



x = tfj + e l9 ^ 2 + e i(t> (* 3 + e*tf 4 ) 



(3.19) 



where 9,(j),£ are three U(l) phases. They parameterize the various boundary conditions 
compatible with self-adjointness of T> w . 
Written explicitly ( |3.19| ) reads, 

/ ^(0) \ / cos(f -f)+e ie cos(f + f) \ 
±«(0) -i(sin(f - I ) + e ie sin(f + f )) 

e**(cos(§ - \ ) + e* cos(f + f )) 



^3.20) 



-ie < *(sin(f - f ) + e* sin(f + f )) / 



Note that the boundary values of ip at the point a = gets contribution from the 
eigenvectors ^1 and ^2 while the boundary conditions at o = R gets contribution from 
and \&4. 

The eigenvalues of the operator T> w are then found by solving by the equation 



a 



— 2ia>tana <9 CT ] ip(a) = \ip(cr) 



(3.21) 



subject to the boundary conditions Q3.26 . 



, p.27p . We seek solutions of the form 



tlj(a) ~ e ipa (3.22) 
Correspondingly, one finds the general solution is a given by 

ip(a) = Ae ip+a + Bj p ~ a (3.23) 

with 

p± = a; tan a ± \ u 2 sec 2 a H — (3.24) 
V a 



Using ( p.24|) in the equation ( |3.20|) we are led to a set of equations which are underde- 



termined. This situation can be remedied by making the choice 8 = £ which allows us to 
determine the eigenvalues A readily : 

\ n = (a--)u 2 + a(— ) , (3.25) 
a V R J 

where n is an integer. 

Note that when 9 — £, the boundary conditions lead to the phase relation 

^(0) = e**i/)(R) 0^(0) = e^d^R) (3.26) 

It is therefore very suggestive to interpret the angle as a Aharonov-Bohm phase factor 
between the end points of the strings, where the quarks are placed. With the hindsights 
obtained from QCD, one can suggest that this phase arises due to chromoelectric field 
stretched between the quarks at the endpoints. However as it is impossible to calculate this 
phase within this model this can only arise in terms of as an self-adjoint extension parameter 
entering at the boundaries. 

Also within this restricted set of boundary conditions one finds also that at the end 
themselves: 

(cos a + e td (l + sin a>)d a ijA + iu (e ie cosa + (1 + sina)) ^ : = (3.27) 
Mixed boundary conditions of this type has appeared in the string literature in the discussion 



of Neumann-Dirichlet tadpoles [0 . If ^ is to be interpreted as a single particle wavefunction 
then one can relate the ratio - to the transverse velocity of the string ends with respect to 
the string center of mass. 



Carrying out the functional sum over all the eigenvalues, following [§] one is immediately 
led to the following potential 

V{R) = 2aR - {D - * )n -=L (3.28) 

. h - (6)2 



Note that the potential again depends on the ratio -. When b is set equal to zero we 
recover the original result due to Liischer et. al. [J7[]. Note that if we hold detr fixed then 
one can think of the tuple (a, b) as a 2d vector on a minkowskian worldsheet whose time 
component is a and the spatial component is ( directed along the length of the string) b. 
Thus the "spatial velocity" would be given by the ratio -. The quantum correction to the 
classical potential then has the correct "Lorentz transformation" factor. Note also that the 
model doesn't require any condition on spacetime dimension D as long as D > 2, unlike the 
treatment in || where one requires D — > oo. This features makes the Siegel's string a more 
general model for hadronic strings compared to the Nambu-Goto string. 

However, from equation ( |3.28| ) it appears that the parameters a, b are arbitrary as long 
as the trace and determinant of rremains positive, i.e. a > and a 2 > b 2 . To stabilize the 
value of - one can introduce potentials for the matrix r mn at the tree level. For instance, if 
we restrict ourselves to quadratic and higher order potentials we can write 

V(r) = 9l Tr(r 2 ) + g 2 det(r) + ■■■ (3.29) 

where gi, g 2 are small coupling constants and the ellipsis represents higher order polynomials 
of r. A short calculation shows that a minima of the potential is obtained for fixed R if g\ 
is negative ( provided a > 0, so that the string tension is positive ). 

In presence of such potentials for r the partonic propagators are modified in the infrared 
regime, specifically picking up higher order poles at p = 0. For instance, if we use the 
quadratic potential in (|3.29| ), the Fock-Schwinger-Dewitt propertime parameterization lead 
to the following propagator ( symbolically) 

J dre-^ p2+9r2+ - ) ~ J dre- Tp \l + gr 2 + ■ ■ •) (3.30) 



1 A 

p 2 2p 



2 ,,,+•■■ (3-31) 



where g represents both the coupling constants in fl3.29| ). 



The appearance of the ^ is interesting from the fact that such singular propagators lead 
to confining potentials. Singular propagators of this type also appear in the dual super- 
conductor picture for QCD [[14] described by Higgs-type models which curiously supports 



stringy Nielsen-Olesen vortices[|T5|. This suggests that if one uses propagators with higher 
order infrared singularities for the partons, the corresponding string theory might contain 
potential terms for the "moduli" field r. 

IV. COMMENTS 

We have examined the Siegel's action for QCD-like strings as an viable model by calcu- 
lating the quark-antiquark potential in this model which is more general than the results of 
Luscher et. al. As stated earlier, in this model the matrix field r mn suffers from the same 
type of runaway problem as the moduli fields in fundamental strings. However, unlike mod- 
uli fields whose potential are prohibited by general covariance or supersymmetry one can 
stabilize the field r mn by introducing a potential for this model. We argued that the role of 
such potentials is to modify the partonic propagators at low-energies. However, such behav- 
ior could be welcomed from the view of dual-superconductor picture for QCD which sustains 
composite stringy solutions [O] like the Abrikosov flux tubes for Type II superconductors. 
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